Finite size scaling in asymmetric systems of percolating sticks 
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We investigate finite size scaling in percolating widthless stick systems with variable aspect ratios in an 
extensive Monte Carlo simulation study. A generalized scaling function is introduced to describe the scaling 
behavior of the percolation distribution moments and probability at the percolation threshold. We show that 
the prefactors in the generalized scaling function depend on the system aspect ratio and exhibit features that are 
generic to whole class of the percolating systems. In particular, we demonstrate existence of characteristic aspect 
ratio for which percolation probability at the threshold is scale invariant and definite parity of the prefactors in 
generalized scaling function for the first two percolation probability moments. 
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I. INTRODUCTION 

Recently there is an increasing interest in the randomly dis- 
tributed stick (rodlike) particles |ll|-|5|], due to promising de- 
velopments in the area of the conducting rodlike nanoparti- 
cle networks, such as carbon nanotubes, silicon, copper, and 
silver nanowires, with applications in electronics [5-7], opto- 
electronics JH], and sensors Jslliol]. On the theoretical side, 
most of the work done until now in the field of percolation of 
random networks is done for lattice percolation fTlUTsll . The 
random stick networks are an important representative of con- 
tinuum percolation [19-22]. Random stick percolation and 
lattice percolation fall into the same universality class having 
the same critical exponents ||2lll . Previous studies established 
that all systems fall on the same scaling function if dimension- 
ality of the system, percolation rule, boundary conditions, and 
aspect ratio are fixed B15I1 . In applications, aspect ratio of the 
rectangular system is the only variable parameter, e.g., geom- 
etry of transistor gate channel in the carbon nanotube transis- 
tors ||4j,|5|]. The objective of the present paper is to describe 
in consistent way finite size scaling of average percolation 
density and standard deviation for the asymmetric rectangu- 
lar stick systems with free boundaries. From general scaling 
arguments one would expect that for all finite size systems 
their convergence is governed by an exponent —1/v iflal . For 
two-dimensional (2D) systems v = 4/3 fl8ll . Following ini- 
tial Ziff's publication [13], Hovi and Aharony HE] argued 
that the irrelevant scaling variables in the renormalization- 
group treatment of percolation imply a slower leading-order 
convergence of percolation probability to its infinite-system 
value, characterized by an exponent —1/v — 8, whose value 
was deducted from the Monte Carlo work of Stauffer to be 
6 » 0.85 HI. Further it was showed that for lattice perco- 
lation on the square system leading exponent of the average 
concentration at which percolation first occurs is —1/v — 6, 
where 9 » 0.9 01711 . All the previous studies were performed 
for symmetric systems. We show that only in symmetric case 
the exponent of average density is — 1/v — 6. In asymmetric 
systems, we observe a leading —1/v exponent. Another quan- 
tity, the percolation probability at the percolation threshold in 
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symmetric bond percolating systems is size independent, i.e., 
scale invariant 112311 . Until now, this behavior is not observed 
in other types of random percolating systems. We will demon- 
strate that asymmetric systems can exhibit scale invariant be- 
havior. 

In this paper, we investigate finite size scaling of the asym- 
metric rectangular stick systems with free boundaries. Both 
from renormalization group considerations and in the simu- 
lations, we find that aspect ratio strongly influences scaling 
behavior of the percolation probability distribution function 
moments, i.e., average density of sticks at which percolation 
first occurs and variance of percolation probability distribu- 
tion function. A generalized scaling function is introduced, 
with aspect ratio dependent prefactors and constant exponents 
of the expansion. Finally, it is shown that percolation prob- 
ability of the asymmetric infinite sticks system at the critical 
threshold density agrees with Cardy's analytic formula 11211 . 



II. NUMERICAL METHOD FOR CALCULATION OF 
PERCOLATION PROBABILITY 

Monte Carlo simulations, coupled with an efficient clus- 
ter analysis algorithm and implemented o n gr i d p latform, are 
used to investigate the stick percolation II24h27I1 . We con- 
sider two dimensional (2D) systems with isotropically placed 
widthless sticks. The sticks of unit length are randomly po- 
sitioned and oriented inside the rectangular system of width 
L x and height L y . Two sticks lie in the same cluster if they 
intersect. System percolates if two opposite boundaries are 
connected with the same cluster. The aspect ratio r is defined 
as the length of the rectangular system in percolating direction 
divided with the length in perpendicular direction. We define 
the normalized system size as a square root of the rectangu- 
lar area L = <J L x L y (geometric average), which represent 
the length of the square system with the same area. The be- 
havior of stick percolation is studied in terms of the number 
stick density n = N/ L 2 . The percolation threshold of infinite 
system is defined by the critical density n c » 5.63726 l27ll . 
Monte Carlo simulations are performed for wide range of 
the aspect ratios, 0.1 < r < 10. In order to ensure the 
same precision for small and large systems we collected more 
than Nmc — 10 9 Monte Carlo realizations for small sys- 
tems L < 10, down to Nmc — 10 7 for the largest system 
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L = 320. Using appropriate functions for fitting data and 
the least squares fitting methods excellent fits were obtained 
(R 2 > 0.9999) for all analyzed systems with L > 16. Statis- 
tical errors for the calculations are estimated in conventional 
fashion using standard deviation l24ll . 

Percolation probability function Rn.l.t, is the probability 
that the system with N sticks, size L and aspect ratio r per- 
colates. It is convenient to pass from the discrete percolation 
probability function Rn.l.t for TV sticks to a probability func- 
tion for arbitrary stick density n l27ll 



m even. Expanding the percolation probability function near 
the critical point we find that 



F{x,{y i },z) = F(0,{0},0) + f Q (x,z)+^f i {x,z)y i + ... 



i=i 



(6) 



where the functions fo(x, z) and fi(x, z) are defined by 
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with percolation probability distribution function defined as 
Pn,L,r = dRn.L.r/dn. The average stick percolation den- 
sity at which, for the first time, a percolating cluster connects 
boundaries of the system is 

(n)L,r = / nP n ,L,rdn = — ^ (1 - R N x,r), (2) 
•*° N=a 

where the last equality follows from integrating by parts. An- 
other important parameter of probability distribution function, 
P n ,L,r, is variance A| r = (n 2 ) Lyr - (n)| r , where (n 2 ) L , r 
is calculated as 



3,1=0 



dx^dyidz 1 



SP z , for j + I even. (8) 



Since percolation probability distribution function P n .L,r = 
dR n ,L,r/dn gives the probability distribution for a system of 
size L and aspect ratio r to percolate for the first time at stick 
density n, we can define the moments of this distribution 
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/ n 2 P n:L:r dn= — y2(N+l) (1-Rn.l.t)- (3) Substituting Eqs. ©-® in Eq. © the fcth moment scales as 



Equations (f2]i and (f3]l allow calculations of the first two mo- 
ments directly from discrete percolation probability function 
Rn,l,t- This is computationally more efficient since it avoids 
calculation of function R n .L.r with high resolution. 



III. GENERALIZED SCALING FUNCTIONS FOR 
MOMENTS 

The percolation probability function is related to the uni- 
versal scaling function B15I1 



^k({Vi},z) = L 



- Tr k l v 



(10) 



where we introduce general functions g. For odd k, go(z) 
is odd function and gi (z) are even functions of z. For even k, 
go(z) is even and gi (z) are odd functions. Therefore, observed 
parity of prefactors in respect to z should be independent of 
the type of the system. 

From Eq. ( fTOb the scaling behavior of the (n)L.r can be de- 
scribed with generalized moment scaling function with aspect 
ratio dependent coefficients 



Rn,L,r = F(x, {j)t},z). 



(4) 



The arguments of the universal scaling function F are x — 
A(n — n c )L 1 / 1/ , jji = BiLUiL~ di , and z = Cln(r), where 
A, {Bi}, and C are the nonuniversal metric factors, {wj} are 
the irrelevant variables, and {Qi} are the corrections to scal- 
ing exponents, (i = 1,2,...). Using free boundary conditions 
and considering two complementary systems: the sticks and 
empty space around the sticks, we can conclude that either the 
sticks percolate in one direction or the empty space percolates 
in the opposite direction: 



{n) L ^ = n c + L-^ v Y d a i {r)L- 



(11) 



where {9i} are the corrections to scaling exponents. The 
zeroth-order correction to exponent 9q should be zero 111 811 . In 
analogy to {ri)L,r, for variance A| r we introduce expansion 



^ r = L- 2 /»Y;bdr)L-°>. 



(12) 



F{x, {fji}, z) + F(-x, {-&}, -z) = 1. 



(5) 



Taking the derivative with respect to x, jji, and z and evaluat- 
ing the derivatives at x = iji = z = 0, (i.e., at n = n c , L — > 
oo, r = 1), we conclude that d m F/dx 3 dyi ... dz l I = 0, for 



From Eq.[l0]and the parity of go(z) and gi(z), for the zeroth- 
order and the first-order prefactors for (n)L,r an d ^1 r near 
ln(r) = (i.e., z = 0), we obtain approximate expressions for 

Oo(r) ~ a , ln(r) + a ,i ln 3 (r), ai(r) w a lf0 + cti.i ln 2 (r), 
bo(r) w b . + b lfi ln 2 (r), 6i(r) « b lfl ln(r) + bi A ln 3 (r). 
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IV. RESULTS AND DISCUSSION 

The results for percolation probability R n ,L,r and distribu- 
tion P nt L r function are shown in Fig.Q] One observes that 
the slope of percolation probability function increases with the 
increase of the system size. The percolation probability func- 
tion curves intersect approximately at n c . The fine behavior of 
percolation probability at n c will be discussed bellow. With 
the increasing system size, the standard deviation of proba- 
bility distribution function decreases to zero. Also, average 
stick percolation density (n) L,r, which corresponds roughly 
to maximum of probability distribution function P n .L.r, ap- 
proaches to percolation threshold n c . For r < 1, (n)L,r 
converges to n c from below with increase of the system size 
L. The reason for this is that narrow finite systems will be 
spanned already at lower densities than n c . For r > 1, (n)i, r 
converges from above, while for symmetric systems (r = 1) 
is roughly centered at n c , see FigQ] From Fig. [2] one can see 
that average stick percolation density (n)L,r for aspect ratio 
higher than one is monotonically decreasing function of the 
system size L. Somewhat surprising, for aspect ratios lower 
than one, {n)L >r is not monotonic function and has local min- 
imum, i.e., for small systems (n)L,r is a decreasing function, 
which passes through n c , reaches a minimum and after that 
converges to n c from below. In inset of Fig. [2] one can see 
that for large system sizes all the curves show power law con- 
vergence to the percolation threshold n c with exponent —1/v, 
except in the symmetric case, i.e., r = 1, where exponent is 
—1/v — 6\. Absolute values of the leading-order prefactors 
are the same for aspect ratios r and 1/r. Higher exponent of 
symmetric systems comes from the basic physics of percola- 
tion that is connectedness. We can illustrate this using a sim- 
plified image of site percolation by introducing quantity R(p) 
as probability that the sites with occupancy p form a span- 
ning path. The percolation probability R(p) and occupancy p 
are equivalent to the percolation probability function R n ,L,r 
and stick density n, respectively. In this image, a cell coming 
out of the renormalization (coarse graining) transformation is 
occupied only if it contains a set of sites that span this cell. 
The universal scaling function in previous section reflects the 
fact that the probability of the spanning system at the percola- 
tion threshold R(p c ) remains unaltered under this transforma- 
tion 01811 . Therefore fixed point of this system, i.e., the critical 
percolation threshold, p c is satisfying relation p c = R(p c )- 
We can expand the percolation probability around the perco- 
lation threshold p c , \R(p) — R(p c )\ dR/dp\ Pc \p—p c \. Also, 
if we renormalize lattice by a length factor b, close to p c , char- 
acteristic length changes as £/&. Since £ ~ \p — p c \ ~ 7 , we 
can write another relation, \R(p) — R(p c )\~' f « \p — p c \ ' /b, 
connecting characteristic lengths before and after renormal- 
ization. From these two relations one can conclude that criti- 
cal exponent should be 



-l/ 7 



In dR/dp\ Pc 
In 1/6 



(13) 



From Fig[TJ one can see that probability density P nc ,L,r which 
is derivative of R n x,r at n c , is always larger for symmetric 
systems than for asymmetric systems of the same size. There- 



fore, from Eq. ( fT3l . one expects higher absolute value of the 
exponent in symmetric compared to asymmetric systems. An- 
other conclusion one can draw from this analysis, is that the 
observed exponents are a result of the interplay of the charac- 
teristic length and the system shape. Usually, such behavior 
is attributed to a competition between two-dimensional and 
three-dimensional (or one-dimensional and two dimensional), 
e.g., in Ising model for slab geometries, cf. Ref B28I1 . In this 
system we observe that there is sharp transition in nature of 
scaling when we pass from symmetric to asymmetric system, 
and a competition between exponents characteristic for sym- 
metric and asymmetric systems. 
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FIG. 1. (Color online) Percolation probability function R n ,L,r (a, b, 
c) and probability distribution function P n ,L,r (d, e, f) for stick per- 
colation on rectangular systems with free boundary conditions and 
increasing system size from L — 20 to 200 for three aspect ratios 
r — 0.5, 1, and 2. The direction of the increase of L is indicated on 
graphs. The vertical dashed lines denote the value for the percola- 
tion threshold n c , while the horizontal dashed lines (a, b, c) denote 
the percolation probability on the infinite systems at the threshold 

Rn c ,L — >oo,r • 

From Monte Carlo simulation data we have obtained the 
two leading-order terms of (n)L,r in Eq. (fTTT i by interpolation, 
cf. Ref. 12911 . Results of the analysis are shown in Fig. [3] The 
zeroth-order prefactor is zero for symmetric system r = 1, 
and it is odd function on a logarithmic scale, i.e., ao(r) — 
—ao(l/r). We have verified the obtained results by interpola- 
tion through symmetrizing data points ((n)z,,r+ (n)i ! i/ r )/2. 
The fitting coefficients for prefactors a^j are calculated us- 
ing the least squares fitting methods and given in Table U cf. 
Fig. [3] Influence of higher order terms were comparable or 
smaller than simulation data error and we could not extract 
them with sufficient precision. For the first-order correction, 
we obtain 6 X = 0.83(2)forr = l,cf. Ref. OS. The resid- 
ual aspect ratio dependence of 6\ cannot be further analyzed 
without provision of retaining the first two terms in Eq. ( fTTT ). 
The system size where the average density reaches minimum 

isi min w (-ai(r)/ao(r)(l + i/0i)) 1/e \ cf. Eq. O- For 
narrow systems, r < 1, L m \ n diverges as l/ln(r) as r ap- 
proaches one. For L < L m ; n the first-order term is dominant. 
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FIG. 2. (Color online) The dependence of the average stick per- 
colation density (n)_t. r on the system size L and aspect ratio r. 
The values are obtained from Monte Carlo simulations and calcu- 
lated using Eq. {2). The values are given for aspect ratios r — 
0.7, 0.8, 0.9, 0.95, 0.98, 1 (solid lines) and their inverse values r = 
1/0.7, 1/0.8, 1/0.9, 1/0.95, 1/0.98 (dashed lines). The horizontal 
bold line denotes the expected value for the percolation threshold n c . 
Inset: The same data is shown in logarithmic scale to demonstrate the 
same power law convergence of the r and 1/r pairs. 
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FIG. 3. (Color online) Prefactors (a) and exponents (b) are shown for 
the two leading-order terms of generalized scaling function for aver- 
age stick percolation density (n) l,t, Eq. dl lb . The zeroth-order pref- 
actor is odd function on a logarithmic scale, i.e., ao(r) = — aoil/r), 
and the zeroth-order exponent is —1/v (solid lines). The first-order 
prefactor is even function, i.e., ai(r) = txi(l/r) and the first-order 
correction to the scaling exponent is 8\ — 0.83(2) for r = 1 (dashed 
lines). 



The variance prefactors and exponents for the two leading- 
order terms are shown in Fig. |U The prefactors and expo- 
nents are obtained by fitting A 2 r with the first two terms in 
Eq. ( [TZi i. cf. Ref. IBOll . The fitting coefficients are given 
in Table U and obtained prefactor dependences on r are given 
in Fig. [4] The zeroth-order prefactor of A| is even function 
on a logarithmic scale, i.e., &o(r) = bo(l/r), as one can see 



TABLE I. Results for the coefficients dij and bi t j, where i,j £ 
{0, 1}. The results are obtained using the least squares method. 
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from a coarse observation of the percolation probability dis- 
tribution function in Fig. Q] Asymmetry of the variance, i.e. 
A| r ^ A 2 L , /r , is the first-order effect, cf. Eq. ( fT2b . 
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FIG. 4. (Color online) Prefactors (a) and exponents (b) are shown 
for the two leading-order terms of generalized scaling function for 
variance Af, r , Eq. (121 l. The zeroth-order prefactor is even function 
on a logarithmic scale, i.e., bo(r) = 6o(l/r) and the zeroth-order 
exponent is —1/v (solid lines). The first-order prefactor is odd func- 
tion, i.e., b\(r) — — 6i(l/r), and the first-order correction to the 
scaling exponent is 8i = 0.80(5) for r = 1 (dashed lines). Prefac- 
tor &i(r) passes through zero for r between 1/0.9 and 1/0.8 causing 
higher error bars of 8i . 

Finally, we investigate scaling behavior of percolation 
probability at the percolation threshold (Rn c ,L,r)> using 
the generalized scaling function R nc ,L,r = Rn c ,L^oo,r + 
ci(r)/i + c 2 (r)/L 2 , cf. Ref. JH. "The results for prefac- 
tors ci(r) and c 2 (r) are shown in Fig. [5] For two limiting 
cases (r < 0.1 and r > 10), the prefactors are close to zero, 
what is consistent with behavior observed in Fig.Q] Between 
these two limiting cases, one can observe that both prefactors 
are close to zero for r = 2.25(5). Furthermore, at this as- 
pect ratio, we could not observe existence of the higher order 
terms. This means that at the percolation threshold percola- 
tion probability is the same for all the systems, independent 
of system size L, having the value R nc .ra2.25 ~ 0.135. Scale 
invariance, i.e., R nc ,L,r = Rn c ,oo,r is already seen and intu- 
itively understood for bond percolation in symmetric systems, 
where R Vr .=o.5, r =i = 1/2 is independent of system size L, 
cf. Ref. 12311 . The reason for observed system size invariance 
of percolation probability at the threshold in the asymmetric 
stick system is existence of multiple zero of at least second 
order at this point in universal scaling function. Regarding 
value of percolation probability at the percolation threshold 
of infinite system, we find that Cardy's analytical model de- 
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rived for the lattice percolation describes simulation data, cf. 
Fig. [5] The deviation between analytical value for lattice and 
Monte Carlo value for stick percolation is less than statisti- 
cal error of the simulation data, i.e., less than 10~ 5 . Per- 
colation probability for 2D sticks system therefore satisfies 
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FIG. 5. (Color online) (a) Prefactors for finite size scaling of the 
percolation probability at the percolation threshold are shown for the 
two leading-order terms, (b) Percolation probability at the threshold 
for infinitely large system -Rn c ,L->oo,r- Points represent Monte Carlo 
data for stick percolation, while line represents results of Cardy's 
model for the lattice percolation. The error bars are much smaller 
than the size of the symbols for all r 



presented methodology could be used to model accurately 
these properties for any percolating system. We find that 
the zeroth-order prefactor of average percolation density is 
odd function with respect to ln(r). This explains the faster 
convergence of average percolation density for symmetric 
systems than expected from general scaling arguments. We 
also observe that there is a characteristic aspect ratio for 
which percolation probability at the percolation threshold is 
system size independent. In addition, for infinite system, we 
find that percolation probability at the critical threshold den- 
sity shows excellent agreement with the Cardy's prediction 
for lattice percolation. The presented results confirm that 
continuum percolation belongs to the same universality class 
as lattice percolation in the sense that value of percolation 
probability at the threshold for infinitely large systems is the 
same for lattice and continuum percolation. One should note 
that a number of other features observed in this work should 
be a common characteristic within the class, e.g., existence 
of the aspect ratio where the percolation probability at the 
threshold is scale invariant and parity of the moments of the 
percolation probability distribution function. This opens up 
the question of the particle shape influence on prefactors, 
weather it is possible to find systems where the observed 
behaviors are more pronounced, and finally the question of 
the general form of the prefactors for describing different 
systems. 
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[29] In order to obtain 9o with high precision successive slopes be- 
tween points L and L/2 11711 are calculated. The slope is de- 
fined as (ln|(n)_L jr — n c | — \n\{n) L / 2 ,r ~ w c |)/ln2. This 
is approximately (ai (r)/a (r)/ In 2)(1 - 2 Bl - e o)£-ei+fo _ 
1/v — 6o when |ao(r)| S> \ai(r)\L~ 9l+0 ° . Therefore, succes- 
sive slopes should fall on a straight line when plotted as a func- 
tion of L~ 9l+ with intercept at —1/v — 8q. The intercept is 
not highly sensitive to the value of — 9i + 9o- Note that expo- 
nent 6a can not be calculated using successive slopes if ao (r) 
is close to zero. 
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|&i(r)|I/ _9l+fl ° is always satisfied. The results for the 
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symmetrizing data points (A|, jT . — A|_ 1 / r )/2. 



